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We examine the nonlinear dynamics of a closed Friedmann-Robertson- Walker universe in the 
framework of Brane World formalism with a timelike extra dimension. In this scenario, the Fried- 
mann equations contain additional terms arising from the bulk-brane interaction which provide a 
concrete model for nonsingular bounces in the early phase of the Universe. We construct a nonsin- 
gular cosmological scenario sourced with dust, radiation and a cosmological constant. The structure 
of the phase space shows a nonsingular orbit with two accelerated phases, separated by a smooth 
transition corresponding to a decelerated expansion. Given observational parameters we connect 
such phases to a primordial accelerated phase, a soft transition to Friedmann (where the classical 
regime is valid), and a graceful exit to a de Sitter accelerated phase. 

PACS numbers: 98.80.Cq, 04.60.Ds 
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I. INTRODUCTION 

Although General Relativity is the most successful the- 
ory that presently describes gravitation, it presents some 
intrinsic crucial problems when we try to construct a cos- 
mological model in accordance with observational data. 
In cosmology, the ACDM model gives us important pre- 
dictions about the evolution of the universe and about 
its current state[l|. However, let us assume that the ini- 
tial conditions of our universe were fixed when the early 
universe emerged from the semi-Planckian regime and 
started its classical expansion. Evolving back such ini- 
tial conditions using the Einstein field equations, we see 
that our universe is driven towards an initial singularity 
where the classical regime is no longer valid 0]. 

Notwithstanding the cosmic censorship conjecture^, 
there is no doubt that General Relativity must be prop- 
erly corrected or even replaced by a completely new the- 
ory, let us say a quantum theory of gravity. This de- 
mand is in order to solve the issue of the presence of 
the initial singulariy predicted by classical General Rel- 
ativity, either in the formation of a black hole or in the 
beginning of the universe. While a full quantum gravity 
theory remains presently an elusive theoretical problem, 
quantum gravity corrections near singularities formed by 
gravitational collapse have been the object of much re- 
cent research, from quantum cosmology^ HJ to D-brane 
theory In the latter scenario, extra dimensions are 
introduced constituting the bulk space. In the case of 
spatially homogeneous and isotropic cosmologies, the ba- 
sic resulting distinction between the two approaches lies 
in the corrections introduced in the Friedmann Hamil- 
tonian constraint, leading either to modifications in the 
kinetic energy terms or to extra potential energy terms. 
In both cases we may have bounces in the scale factor 
corresponding to the avoidance of a singularity in the 
models. In this context, the initial conditions from which 



our universe has evolved should depend crucially on the 
adopted version of the theory to describe the dynamics 
around the singularity. 

One of the most important characteristics of our uni- 
verse, supported by observational data, is its large scale 
of homogeneity and isotropy. In fact, the scale of ho- 
mogeneity and isotropy is empirically well accepted for 
distances above 100 Mpc. Indeed this is the main rea- 
son that makes the geometry of Friedmann-Lematre- 
Robertson- Walker (FLRW) a powerful theoretical tool 
for the construction of a cosmological scenario pj. How- 
ever, when we consider a homogeneous and isotropic 
model filled with baryonic matter, we find several dif- 
ficulties when we take into account the primordial state 
of our universe. Among such difficulties, we can mention 
the horizon and flatness problems [lj. 

As a possible solution to theseproblems emerged the 
so called Inflationary Paradigm [l], [l(| ■ Although this fun- 
damental paradigm allows to solve the horizon and flat- 
ness problems, inflationary cosmology does not solve the 
problem of the initial singularity. Therefore, nonsingular 
models from a new theory which provide alternative so- 
lutions to these problems should be strongly considered. 

In this paper we adhere to the brane world scenario, 
where a timclikc noncompact extra dimension is intro- 
duced, constituting the bulk space, and all the mat- 
ter content of our universe would be trapped on a 4- 
dimensional spacetime embedded in the bulk. At low 
energies General Relativity is recovered Q, but at high 
energy scales significant changes arc introduced into the 
gravitational dynamics and the singularities can be re- 
moved Q. 

While spacelike extra dimensions theories have re- 
ceived more attention in the last decades @, stud- 
ies involving extra timelike dimensions have been 
considered ll| despite from the fact that propagating 
tachyonic modes or negative norm states may arise 
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due to timelike extra dimensions. These modes have 
been regarded as problematic once they might violate 
causality fl2l| by considering interactions among usual 
particles. Issues like the exceedingly small lower bound 
on the size of timelike extra dimensions (l3j . the imagi- 
nary self-energy of charged fermions induced by tachyonic 
modes - which seems to cause disappearance of fermions 
into nothing - and the spontaneous decay of stable parti- 
cles induced by tachyons with negative energy are major 
difficulties [l2|]. Nevertheless, in order to address the cos- 
mological constant problem in Kaluza-Klein theories flij 
or reconcile a solution of the hierarchy pro blem with the 
cosmological expansion of the universe [15|, timelike ex- 
tra dimensions have been considered. On the other hand, 
it has been shown in [l(| that the appearance of mass- 
less ghosts in an effective four-dimensional theory can 
be avoided by considering topological criteria in Kaluza- 
Klein theories with extra compactified time-like dimen- 
sions. Moreover, avoidance of propagating tachyonic or 
negative norm states can also be achieved by considering 
a noncompact timelike extra dimension^, which is the 
case in the model of this paper. 

We organize the paper as follows. In section II we 
present a brief review of the modified Einstein field equa- 
tions in the Branc World scenario. In section III, we con- 
struct a nonsingular cosmological scenario sourced with 
dust, radiation and a cosmological constant. In section 
IV, we show that given the observational parameters, we 
can connect such phases to a primordial accelerated ex- 
pansion, a soft transition to Friedmann (where the clas- 
sical regime is valid), and a graceful exit to a de Sitter 
accelerated phase. As our final remarks, we discuss some 
of its possible imprints in the physics of cosmological per- 
turbations. 
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In the above ^R is the Ricci scalar of the Lorentzian 
5-dim metric qab on A4, and ^'R is the scalar curva- 
ture of the induced metric g a b on E. The parameter a is 
denoted the brane tension. The unit vector n A normal 
to the boundary S has norm e. If e = — 1 the signature 
of the bulk space is (—,—,+,+,+), so that the extra 
dimension is timelike. The quantity K = K a b g ab is 
the trace of the symmetric tensor of extrinsic curvature 
K a b = Ya C y,b D VcriB, where Y A (x a ) are the embed- 
ding functions of S in .A/fjlT]]- While Lj(q a b-, p) is the 
Lagrangean density of the perfect fluid [la] (with equa- 
tion of state p = ap), whose dynamics is restricted to 
the brane E, L5 denotes the lagrangian of matter in the 
bulk. All integrations over the bulk and the brane are 
taken with the natural volume elements \J —e ^g d 5 x 
and \J — ( 4 )</ d 4 x respectively. K5 and K4 are Einstein 
constants in five and four-dimensions. We use units such 
that c = 1. 

Variations that leave the induced metric on S intact 
result in the equations 
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while considering arbitrary variations of gAB an d taking 
into account © we obtain 
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II. FIELD EQUATIONS 



For the sake of completeness let us give a brief intro- 
duction to braneworld theory, making explicit the spe- 
cific assumptions used in obtaining the dynamics of the 
model. We rely on references [H, , and our notation ba- 
sically follows @]. Let us start with a 4-dim Lorentzian 
brane £ with metric g a b, embedded in a 5-dim confor- 
mally flat bulk M with metric gAB ■ Capital Latin indices 
range from to 4, small Latin indices range from to 3. 
We regard £ as a common boundary of two pieces M.\ 
and M.2 of M and the metric g a b induced on the brane 
by the metric of the two pieces should coincide although 
the extrinsic curvatures of £ in Mi and M.2 are allowed 
to be different. The action for the theory has the general 



where S a b = K a b — Kg a b, and r a b is the energy mo- 
mentum tensor associated to L4. In the limit K4 — >■ 
00 equation ([3]) reduces to the Israel-Darmois junction 
condition fijl 
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We impose the ^-symmetryQ and use the junction con- 
ditions ([4]) to determine the extrinsic curvature on the 
branc, 
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Now using Gauss equation 
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together with equations ([2]) and ([5]) we arrive at the in- 
duced field equations on the brane 
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E ab is the electric part of the Weyl tensor in the bulk 
induced in the brane, Tab is the energy-momentum in 
the bulk, and Gn defines the Newton's constant on the 
brane. For a timelike extra dimension we have that 
e = — 1 in our conventions implying that a < in ac- 
cordance with observations. We also remark that the ef- 
fective 4-dim cosmological constant might be set zero, or 
made conveniently small, in the present case of an extra 
timclike dimension by properly fixing the bulk cosmolog- 
ical constant as A5 ~ a 2 . It is important to note 
that for a 4-dim brane embedded in a conformally flat 
empty bulk we have the absence of the Weyl conformal 
tensor projection E ab and of F ab in Eq. ©. 
Accordingly, Codazzi's equations imply that 



resulting in 
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where V a is the covariant derivative with respect to the 
induced metric g ab . Equations ([7]) and (|13|) are the dy- 
namical equations of the gravitational field on the brane. 



III. THE MODEL 

Let us consider a FLRW geometry on the four- 
dimensional brane embedded in a five-dimensional de- 
Sitter bulk with a timelike extra dimension (e = — 1)@. 
Considering comoving coordinates on the brane, the line 
element is given by 
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where a(t) is the scale factor, k is the spatial curvature 
and d£l 2 is the solid angle. 




FIG. 1: Illustration of the behavior of V(a) for parameters 
k = 0.8, A 4 = 1.5 E rad = 0.1, and for E dust = 0.001 (con- 
tinuous line) Edust = 0.09 (dashed line). The increase of the 
dust content for a fixed k excludes the presence of perpetually 
bouncing solutions. 



The matter content of the model, restricted to the 
brane, is given by noninteracting perfect fluids, namely, 
dust and radiation, with respective equations of state 
Pdust = 0, Prad = Prad/3, and energy momentum tensor 



T ab :=T ab ,„.,. t + r ab rad satisfying V b r ab dust = = V h r 
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In this situation we have that 
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+ -(pdust + Prad)(pdust + Prad) ^Qij , (15) 

and Codazzi's equations (fl~2"|) imply that V a IIj = 0, 
consistent with the contracted Bianchi's identities in ([7]) 
and Codazzi's equation (fT3| . The modified Friedmann 
equations have the first integral 



H 



A4 
3 



8ttG 



N 



(Pdust + Prad) 



AirG 



N 



3d 



(P 



dust 



PradY 



(16) 



where H := a/a with d = da/dt. It is worth noting 
that the bounce is solely engendered due to the presence 
of a timclike extra dimension which induces the minus 
sign in the last term of (|16[) . By assuming a spacelike 
extra dimension, we would obtain a plus sign instead that 
provides a singular model. 
Expressing 
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FIG. 2: Phase portrait of the dynamics with the critical point 
Pi (center) and P2 (saddle). Orbits in region II are solutions 
of one-bounce universes with a graceful exit to an accelerated 
(inflationary) phase along the separatrix S. 



where Edust and E ra d are constants of motion, the first 
integral of motion (|16[) can be expressed as the Hamilto- 
nian constraint 



■H = p l+V{a) 



0. 



(18) 



large E ra d the potential V(a) presents no local maximum 
or minimum |20j|. 

The critical points in the phase space are stationary 
solutions of ((20)) . namely, the points of the phase space 
(a = a cr it,p a = 0) corresponding to the zeros of the right- 
hand-side of (f2"0"]) . Here, a crit stands for the real positive 
roots of dV/da. By considering the case of closed geome- 
tries (k > 0), it is not difficult to verify that, depend- 
ing on the values of the parameters (A4, \ a\, E ra d, Edust), 
there are at most two critical points associated with one 
minimum and one local maximum of V(a). In this case, 
the minimum of the potential corresponds to a center 
while the maximum corresponds to a saddle. This con- 
figuration allows us to obtain different types of orbits that 
describes the evolution of universes in this model. In Fig. 
2 we illustrate the phase space portrait of the model for 
A 4 = 1.5, a = 6000, E rad = 0.15 and E dus t = 0.05, and 
for varying k. The value of Edust is sufficiently bounded 
so that V(a) has a well. The critical points P\ (center) 
and P2 (saddle) correspond to stable and unstable Ein- 
stein universes. Typically the model allows for the pres- 
ence of perpetually bouncing universes (periodic orbits 
in region I), and one-bounce universes (region II). Re- 
gion I is bounded by the separatrix S emerging from the 
saddle P 2 . A separatrix also emerge from P 2 towards the 
dcSitter attractor at infinity, defining a graceful exit of 
orbits in region II to an (inflationary) accelerated phase. 
From now on we will restrict ourselves to the case of 
closed geometries. In the next section we will exam what 
kind of orbit would be generated when one considers the 
observational values of (A4, E ra d, Edust)- 
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From (|18[) we derive the dynamical system 
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It is worth noting that the last term in the potential (jT9|) 
acts as an infinite potential barrier and is responsible for 
the avoidance of the singularity a — 0. These potential 
corrections are equivalent to fluids with negative energy 
densities. This is in accordance with the fact that quan- 
tum effects can violate the classical energy conditions, 
and may avoid curvature singularities where classical gen- 



eral relativity breaks down[21|. Such violations tend to 



occur on short scales and/or at high curvatures, which is 
the case of the present models. 

The behavior of the potential V(a) is illustrated in Fig. 
□for k = 0.8, A 4 = 1.5, E rad = 0.1, and for E dus t = 0.001 
(continuous line) Edust = 0.09 (dashed line). We can see 
that the increase of the dust content for a fixed k excludes 
the presence of perpetually bouncing solutions by driving 
the maximum out of the physical space. For a sufficiently 



IV. OBSERVATIONAL COSMOLOGY 

As observational cosmology asserts, the domain of ho- 
mogeneity and isotropy of our present Universe is well 
accepted for scales around the present horizon, which is 
given by ao ~ 10 28 cm (here, the subscript denotes 
the present epoch). In this case, we obtain the following 
observational parameters 



A 4 ~ 1.34 x 10~ 56 cuT 2 
Edust s 2.6 x 10 54 g , 
E rad ^ 4 x 10 78 g cm, 



the Hubble 

-1 



radius is fixed to Hq 



(21) 
(22) 
(23) 

0.77 x 



where 
10" 28 cm- 
Prom Ref. [iH, the brane tension has a lower bound 
which corresponds to |o"| m j n ~ 10 22 g cnr 3 . That is, 
a star with the Chandrasekhar mass will not form an 
event horizon if the brane tension is smaller than |<r| m j n . 
It turns out that this value furnishes us with a curvature 
scale l c = 1/y/Rb = {\Ja/d) ~ 10 34 lp at the bounce 
(where I p is the Planck lenght and Ri, is the Ricci scalar 
at the bounce). In order to guarantee that l c at the 
bounce is not smaller than 10 3 lp, the brane tension must 
be less than 10 85 g cm -3 . Therefore, we have the follow- 
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FIG. 3: (a) The potential V(a) and the phase space (b) in 
the region that encompasses the critical point eti , considering 
the observational values (|21[) and (|24|l . For \a\ ~ |<r| m in, we 
obtain a\ ~ 10 14 cm. Although this primordial accelerated 
phase does not correspond to usual inflation (0.27 being the 
number of e-folds), it is important to remark that as our uni- 
verse has no beginning of time and the cosmological constant 
is small, the particle horizon before the bounce was already 
bigger than the scales of cosmological interest. In (c) we show 
the behavior of the Hubble factor as function of redshift in a 
neighborhood of the bounce given the normalization ao = 1 
and the parameters (f3"0")l . 



FIG. 4: (a)The potential V(a) and the phase space (b) in 
the region that encompasses the critical point ai and com- 
pletes Figs. [3] This critical point is of the order of 10 28 cm, 
coinciding with the domain of homogeneity and isotropy of 
our present Universe, regardless of the value of \a\. In (c) we 
show the behavior of the Hubble factor as function of redshift 
in a neighborhood of the saddle a% given the normalization 
ao — 1 and the parameters (|30l) . The domain > z > — 1 of 
the H > branch corresponds to the final acceleration phase 
approaching de Sitter as z — ¥ — 1, with H — const, at z = — 1. 
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ing physical domain (not spoiling the nucleosynthesis) for 
the brane tension 



10 22 g/ cm £ kl ^ 10 g/ cm 



(24) 



where we have set c = 1. Feeding the Hamiltonian con- 
strain (|18p with k| m j„ and the parameters (|2ip . we ob- 
tain that the spatial curvature is k ~ 0.002 for |<r| > 

k I mill- 

Considering the lower bound limit for \a\, numerical 
calculations show that the potential V(a) has always a 
local minimum at a\ (corresponding to a center) and a 
local maximum at a 2 (corresponding to a saddle). - cf. 
Figs. 3(a) and 4(a). If we increase \a\ by four orders of 
magnitude, we obtain a value of a\ decreased by one order 
of magnitude. On the other hand, the local maximum a-i 
is of the order of 10 28 cm (cf. Fig. 4(a)) for \a\ > [cr| min 
(regardless of the value of |<r|). We exhibit the behavior 
of V(a) and the phase space (a,p a ) trajectory - for the 
parameters (j2Tj) -(24) and \a\ ~ |cr| m j n - in Figs. 3(a) and 
4(a), and Figs. 3(b) and 4(b), respectively. We should 
note that Figs. 3 and Figs. 4 display the same potential 
V{a) and the same universe phase space trajectory in 
distinct ranges of o, complementing each other. 

Given the redshift relation 
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cquation (fll)|) can be rewritten as 
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and Pq t ' — 3Hq /8irG. By fixing the normalization clq = 
1, we obtain the following parameters according to the 
WMAP 7 year results [23] : 







0.26 , O 0rod ~ 10- 5 , n 0A 0.73. (30) 



Odust 



Substituting these parameters in (|26[) . we obtain ilok — 
0.004. In figures 3(c) and 4(c) we show the behavior of 
the Hubble scale factor H with respect to the redshift z. 

It is remarkable that considering the interval of 62 or- 
ders of magnitude of \a\ (cf. (|24"]l). the trajectory in the 
phase space of the above observable universe belongs to 
region II of the phase space (cf. Fig. [5]) corresponding 
to a one-single-bounce orbit. The part of the trajectory 



starting from (a = ai,p ai < 0) is an initial acceleration 
phase that leads the universe through the bounce and 
ends in (a = ai,p ai > 0), when the universe enters in a 
long and smooth decelerated expansion phase. This pri- 
mordial bouncing accelerated phase does not correspond 
to usual inflation, the number of e-folds being 0.27. Note, 
however, that there is no horizon problem in the model. 
In fact, before the bounce, due to its cosmological con- 
stant dominated contraction from the infinity past until 
<Z2, the particle horizon d p is given by 



gL = 



a j 



-da 



10 28 cm, 



(31) 



if a > <i2 . Therefore the particle horizon is already of the 
order of A -1 / 2 , which is constrained by present observa- 
tions to be of the order of the Hubble radius today. Hence 
there is no horizon problem for the scales of cosmologi- 
cal interest. The decelerated expansion Friedmann phase 
ends in the neighborhood of ai with a graceful exit to a 
late de Sitter accelerated phase. 

From © and (f2"T|) . we see that the parameter A 5 of 
the model must be adjusted in a very precise way. In 
fact, A 5 must be very close to 87rGjvk|, which has the 
minimum value 10 5 cm -2 (see Eq. (|2~4")) ) and increases as 
| a | increases, in order to yield the observed value of A4 
given in Eq. ((2Tj) . This is the usual problematic fine- 
tuning of the cosmological constant, of at least 60 orders 
of magnitude as we have seen above, which the present 
model, at least in this first approach, does not solve. It 
turns out that this is similar to an issue contained in the 
Randall and Sundrum model 241. In this scenario the 



brane is embedded in an anti-dc Sitter (4 + 1) spacctime 
and the fine-tuning relation A5 = — k|<7 2 /6 has to be 
satisfied. It was shown in [25| that the Randall-Sundrum 
model is unstable under small deviations from this fine- 
tuning. As a future investigation we will examine if the 
same happens in our model. 



V. CONCLUSIONS 

In the framework of a Brane World formalism with a 
timclike extra dimension, we have obtained a homoge- 
neous and isotropic bouncing model compatible with all 
observations at the background level. It starts with a 
de Sitter contraction from the infinity past, experiences 
a bounce at very small scales, turning to the usual stan- 
dard expanding decelerating phases of radiation and mat- 
ter domination, and a recent transition to an accelerating 
expansion. The bounce itself is caused by the appearance 
of new terms coming from the extra timelike dimension of 
the bulk in the 4-dimensional Friedmann equation, which 
become important at high curavture scales and avoid the 
cosmological singularity, inducing a gravitational repul- 
sion due to the timelike nature of the extra dimension. 
We have two free parameters: the brane tension a and 
the five dimensional cosmological constant A5 . The brane 
tension can assume a wide variety of values, see Eq. 
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but A5 must be highly fine-tuned to the value of a in or- 
der to yield an effective 4-dimensional cosmological con- 
stant compatible with observations. Hence the model 
solves the singularity problem of the standard cosmologi- 
cal model, together with the horizon and flatness puzzles, 
but it does not solve the cosmological constant problem. 

Our next step will be to perturb the model and investi- 
gate the evolution of cosmological perturbations in such 
cosmological background. Indeed, our work in progress 
shows that, if one imposes an unperturbed de Sitter 
bulk, a numerical treatment of linear hydrodynamical 
perturbation in the universe indicates that the bounce 
has the effect of substantially enhance the perturbations, 
nonetheless these perturbations remain bounded with 
Sp/p <C 1 and Sp/p <C 1. However, a general analysis 
of cosmological perturbations in this scenario demands 
also a perturbed de Sitter bulk. In this case the 5-D 
scalar perturbations will induce fluctuations of the Weyl 



tensor projected on the brane which will modify the per- 
turbed field equations. This is a technical and conceptu- 
ally involved problem which will be investigated in future 
publications. 
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